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Abstract. The Epoch of Reionization (EoR) is the epoch in which hydrogen in the Universe reionize after the "Dark Ages". This 
is the second of two major phase transitions that hydrogen in the Universe underwent, the first phase being the recombination 
era in which hydrogen became neutral at redshift « 1100. The EoR, occurs around z ~ 10 and is probably caused by the first 
radiation emitting astrophysical sources, hence it is crucial to our understanding of when and how the Universe "decided" to 
start forming astrophysical objects and how that influenced subsequent structure formation in the Universe. As such, the EoR 
is related to many fundamental questions in cosmology, galaxy formation, quasars and very metal poor stars; all are foremost 
research issues in modern astrophysics. The redshifted 21 cm hyperfine line is widely considered as the most promising probe 
for studying the EoR in detail. In the near future a number of low frequency radio telescopes (LOFAR, MWA, GMRT and SKA) 
will be able to observe the 21 cm radiation arriving from the high redshift Universe. In this paper I present our current picture 
of the ionization process, review the 21 cm line physics and discuss the challenges that the current generation experiments are 
expected to face. Finally, I discuss the potential of SKA in exploring the EoR and the Universe's Dark Ages. 



1. Introduction 

The last couple of decades have witnessed the emergence of 
an overarching paradigm, the ACDM model, that describes 
the formation and evolution of the Universe and its structure. 
The ACDM model accounts very successfully for most of the 
available observational evidence on large scales. According 
to this paradigm at redshift of 1100, about 400,000 years af- 
ter the Big Bang, the Universe's temperature and density de- 
creased enough to allow the ions and electrons to recombine 
into neutral hydrogen and helium (and a very small percent- 
age of heavier elements). Immediately afterwords, photons de- 
coupled from the baryonic material and the Universe became 
transparent leaving a relic radiation, known as the cosmic mi- 
crowave background (CMB) radiation. To date, the CMB ra- 
diation has been observed in many experiments providing one 
of the most compelling evidences for the CDM paradigm (for 
recent results see the WMAP papers, e.g., Spergel et al. 2007; 
Pageetal. 2007; Komatsu et al. 2010). 

The matter-radiation decoupling has ushered the Universe into 
a period of darkness as its temperature dropped below 3000 K 
and steadily decreased with the Universe's expansion. These 
Dark Ages ended about 400 million years later, when the first 
radiation emitting objects (stars, black-holes, etc.) were formed 
and assembled into protogalaxies. The first objects then began 
to produce a spectrum of radiation, especially ionizing ultra- 
violet photons. After a sufficient number of ionizing sources 
had formed, the temperature and the ionized fraction of the gas 
in the Universe increased rapidly and most of the neutral hydro- 
gen eventually ionized. This period, during which the cosmic 
gas went from being almost completely neutral to almost com- 
pletely ionized, is known as the Epoch of Reionization of the 
Universe. 

The most accepted picture on how reionization unfolds is sim- 
ple. The first radiation-emitting objects ionize their immedi- 
ate surroundings, forming ionized bubbles that expand until 



the neutral intergalactic medium (hereafter, IGM) consumes all 
ionizing photons. As the number of objects increases, so do the 
number and size of the bubbles, until they eventually fill the 
whole Universe. However, most of the details of this scenario 
are yet to be clarified. For example: what controls the forma- 
tion of the first objects and how much ionizing radiation do 
they produce? How do the bubbles expand into the intergalac- 
tic medium and what do they ionize first, high-density or low 
density regions? The answer to these questions and many oth- 
ers that arise in the context of studying the EoR touches upon 
many fundamental questions in cosmology, galaxy formation, 
quasars and very metal poor stars; all are foremost research 
issues in modern astrophysics. A substantial theoretical effort 
is currently dedicated to understanding the physical processes 
that trigger this epoch, govern its evolution, and the ramifica- 
tions it had on subsequent structure formation (c.f., Barkana 
& Loeb 2001; Bromm & Larson 2004; Ciardi & Ferrara 
2005; Choudhury &Ferrara 2006; Furlanetto & Briggs 2006). 
However, and despite the pivotal role played by the EoR in cos- 
mic history, observational evidence on it is very scarce and, 
when available, is indirect and model dependent. 

It is generally acknowledged that the 21cm emission line from 
neutral hydrogen at high redshifts is the most promising probe 
for studying the Dark Ages and the EoR in detail. HI fills the 
IGM except in regions surrounding the ionizing radiation of 
the first objects to condense out of the cosmic flow. Computer 
simulations suggest that we may expect an evolving complex 
patch work of neutral (HI) and ionized hydrogen (HII) regions 
(Gnedin & Abel 2001; Ciardi et al. 2001; Ritzerveld et al. 
2003; Susa 2006; Razoumov & Cardall 2005; Nakamoto 
etal. 2001; Whalen & Norman 2006; Rijkhorst et al. 2005; 
Mellema et al. 2006; Zahn et al. 2007; Mesinger & Furlanetto 
2007; Pawlik & Schaye 2008; Thomas et al. 2009). The cur- 
rent constraint strongly suggest that the EoR roughly straddles 
the redshift range of z ~ 20 - 6. 
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In this contribution I will give a short overview of the present 
picture of the Universe's reionization and review the potential 
of current and future instruments -focusing mostly on LOFAR 
and SKA- in detecting it. The paper is organized as follows. 
In section 2 a review of the current observational constraint is 
given. In section 3 the possible reionization sources are dis- 
cussed. Section 4 presents the physics of the 21 cm line with 
the factors that determine the spin and brightness temperature 
(the observable in these experiments) and show a number of re- 
sults obtained from numerical simulations. In § 5 the expected 
foregrounds, the instrument sensitivity and the predicted power 
spectrum measurement and the skewness are discussed. The 
paper ends with conclusions. 

2. Current observational constraints 

Currently, there are a few strong observational constraints on 
the EoR. The CMB temperature and polarization data ob- 
tained by the WMAP satellite allow measurement of the total 
Thomson scattering of the primordial CMB photons off inter- 
vening free electrons produced by the EoR along the line of 
sight (Page et al. 2007; Spergel et al. 2007; Komatsu et al. 
2010). They show that the CMB radiation has, on its way to 
us, only been damped by ~ 9%, indicating that the Universe 
was mostly neutral for 400 million years and then ionized. 
However, the Thomson scattering measurement is an integral 
constraint telling us little about the sources of reionization, its 
duration or how it propagated to fill the whole Universe. 

The second strong constraint comes from specific features in 
the spectra of distant quasars, known as the Lyman-alpha for- 
est. These features, which are due to neutral hydrogen, indicate 
two important facts about reionization: 1) hydrogen in the re- 
cent Universe is highly ionized with neutral fraction of ~ 10 -4 ; 
2) at redshift 6.5, i.e., about 900 million years after the Big 
Bang, the neutral fraction of hydrogen suddenly increases de- 
marcating the end of the reionization process (Fan et al. 2003, 
2006). Despite these data providing strong constraints on the 
ionization state of the Universe at redshifts below 6.5, they say 
very little about the reionization process itself (although see 
Bolton & Haehnelt 2007). 

A third, albeit weaker, constraint comes from measuring the 
temperature of the intergalactic medium at the redshift range 
4-6. The widths of the hydrogen Lyman-o- absorption lines can 
be used to determine the IGM temperature. This property has 
been used by a number of authors who have measured the IGM 
temperature at z» 4 (Schaye et al. 2000; Theuns et al. 2002). 
The temperature of the diffuse IGM at these redshifts depends 
on its reionization history because the thermal timescales are 
long. Theuns et al. (2002) and Hui & Haiman (2003), have 
used these measurement at redshift about 4 to extrapolate back 
in time where they argue that, for a single-phase reionization 
scenarios, hydrogen reionization occurred below redshift z=9. 
Recently, Bolton et al. (2010) have measured the IGM tem- 
perature at even higher redshift (zx 6) and employed similar 
arguments to conclude that reionization happened at redshift 
lower than 9. As said, this constraint is somewhat weaker than 
the previous two as it requires assumptions about the heating 




Fig. 1. This figure shows a sketch of the likely development of the 
EoR. About 500,000 years after the Big Bang (z ~ 1000) hydrogen 
recombined and remained neutral for a few hundred million years (the 
Dark Ages). At a redshift, z ~ 15, the first stars, galaxies and quasars 
began to form, heating and reionizing the hydrogen gas. The neutral 
IGM can be observed with LOFAR through its redshifted 21cm spin- 
flip transition up to redshift 1 1.5. However, many atmospheric, galac- 
tic and extra-galactic contaminants corrupt the 21 cm signal. 

and cooling rates from the redshift at which the temperature is 
measured all the way back to the reionization redshift. 

A whole slew of possible other constraints currently discussed 
in the literature are either very controversial, very weak or, as 
is often the case, both. Most are very interesting and exciting; 
but can only be investigated reliably with a new generation of 
instruments that are expected to come online within a decade. 

To summarize, from the current observational constraints and 
recent numerical and theoretical models of structure formation 
a very simple picture of the Universe's reionization history has 
emerged. A cartoon of the various phases and the objects fea- 
turing in this simple picture is shown in Figure 1 . 

3. Astrophysical Ionization Sources 

The nature of the first ionizing sources is unclear. In general 
two types of astrophysical sources are considered 11 Population 
III stars or miniquasars, both of which represent hypothetical 
but plausible populations of the first objects in the universe, 
and that are significant sources of ionizing photons. 

First stars: The impact of stellar sources have been stud- 
ied by many authors (e.g., Cen 2003; Ciardi et al. 2003; 
Haiman & Holder 2003; Sokasian et al. 2003; Wyithe & 
Loeb 2003). These studies find in general that in order to pro- 
vide enough ionizing flux at or before z = 15, for the usual 
scale-invariant primordial density fluctuation power spectrum, 

a Other ionization sources, e.g., decaying dark matter particles, are 
possible but less likely 
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one needs Population III stars, which provide about 20 times 
more ionizing photons per baryon than Population II (Schaerer 
2002; Bromm et al. 2001), or an IMF that initially is dominated 
by high mass stars. This is in agreement with recent numeri- 
cal simulations of the formation of first stars from primordial 
molecular clouds that suggest that the first metal-free stars were 
predominantly very massive, m, < 100M o (Abel et al. 2000; 
Bromm et al. 2002). 

Mini-quasars: Miniquasars have also been considered as a 
significant ionising source (e.g., Kuhlen & Madau 2005; Oh 
2001; Madau et al. 2004; Ricotti, & Ostriker 2004a,b; Thomas 
& Zaroubi 2008; Zaroubi & Silk 2005; Zaroubi et al. 2007). 
In view of the correlation between central black hole mass and 
spheroid velocity dispersion (Ferrarese 2002; Gebhardt et al. 
2000), this correlation demonstrates that seed black holes must 
have been present before spheroid formation. Theory and ob- 
servations suggests that the seeds from which the super massive 
black holes amounted to at least 1000M o and were in place be- 
fore z ~ 10 (Fan et al. 2003, 2006; Madau & Rees 2001; Silk 
&Rees 1998). 

Decaying Dark Matter particles: In the literature, the influence 
of dark matter particles annihilation or decay on reionization is 
considered. The consensus here is that, given the constraints 
that the CMB impose on the value of Thomson optical depth, 
these process are unlikely to substantially contribute to the 
Universe's ionization. Instead many authors consider whether 
the influence of annihilation/decay of dark matter particles on 
reionization could be used, despite its small amplitude, to con- 
straint the properties of these particles. 

4. The 21 cm Spectral Line as a Probe of the EoR 

In recent years it has become clear that the 21 cm line can be 
used to probe the neutral IGM prior to and during reionization. 
This hyperfine transition line of atomic hydrogen (in the ground 
state) arises due to interactions between and the electron and 
proton spins (Hogan & Rees 1979; Scott & Rees 1990; Madau 
et al. 1997). The excited triplet state is a state in which the 
spins are parallel whereas the spins at the lower (singlet) state 
are antiparallel. The 21 cm line is a forbidden line with a prob- 
ability of 2.9 x 10 _15 i _1 corresponding to a life time of about 
10 7 years. Despite its low decay rate the 21 cm transition line is 
one of the most important astrophysical probes, simply due to 
the vast amounts of hydrogen in the Universe (Ewen, & Purcell 
1951; van de Hulst 1945; Muller & Oort 1951) 

4.1. The 21 cm Spin Temperature 

The intensity of the 21 cm radiation is controlled by one param- 
eter, the so called spin temperature, T spin . This temperature is 
defined through the equation, «i/«o = 3 exp(-T,/T spin ), where 
n\ and no are the number densities of electrons in the triplet and 
singlet states of the hyperfine levels, and T* = 0.068 1 K is the 
temperature corresponding to the 21 cm wavelength. 
In his seminal papers George Field (Field 1958, 1959) used the 
quasi-static approximation to calculate the spin temperature, 
T sp in, as a weighted average of the CMB temperature, Tcmb, 



the gas kinetic temperature, T kin , and the temperature related 
to the existence of ambient Lyman-a photons, T a (Wouthuysen 
1952; Field 1959). For almost all interesting cases one can 
safely assume that 7\ !n = T a (Field 1958; Furlanetto & Briggs 
2006; Madau et al. 1997). The spin temperature is given by: 

T _ T C MB + y'kinTkin + y a Tkin m 
* spin — -, > U ) 

l + ytin + y a 

where yfa„ and y a are the kinetic and Lyman-a coupling terms, 
respectively. The kinetic coupling term is due to collisional ex- 
citations of the 21 cm transitions. The Lyman-a coupling term 
is due to the so called Lyman-c pumping mechanism - also 
known as the Wouthyusen-Field effect - which is produced by 
photo-exciting the hydrogen atoms to their Lyman transitions 
(Field 1958, 1959; Wouthuysen 1952). Finally, it is important 
to note that for the 21 cm radiation to be observed it has to at- 
tain a different temperature than that of the CMB background 
(Field 1958, 1959; Hogan & Rees 1979; Wouthuysen 1952). 

Since the decoupling mechanisms can influence the spin tem- 
perature differently, it is important to explore the decoupling 
issue for various types of ionization sources. For instance, 
stars decouple the spin temperature mainly through radiative 
Lyman-a pumping whereas mini-quasars decouple it through 
a combination of collisional Lyman-a pumping and heating 
(Chuzhoy et al. 2006; Zaroubi et al. 2007), both produced 
by energetic secondary electrons ejected due to the miniqso 
x-ray photons (Shull & van Steenberg 1985). The difference 
between the spin temperature decoupling patterns between the 
two will eventually help disentangle the nature of the first ion- 
ization sources (Thomas & Zaroubi 2008). 

Figure 2 shows the expected global evolution of the spin tem- 
perature as a function of redshift. The blue solid line represents 
Tcmb, which drops as 1 + z. The green line shows the gas tem- 
perature as a function of redshift. At z ~ 200 the gas temper- 
ature is still coupled to the CMB due to Compton scattering 
off residual electrons leftover from the recombination era. At 
redshift ~ 200, however, the gas decouples from the CMB and 
starts adiabatically cooling as a function of (1 + z) 2 until the 
first objects start forming and heating up the gas at redshift be- 
low 30. The spin temperature (shown by the red lines) has a 
somewhat more complicated behavior, at z > 30 it is coupled 
to the gas temperature due to collisional coupling caused by 
a residual electrons leftover from recombination. At z ~ 100 
the efficiency of collisional coupling to the gas drops due to 
the Hubble expansion. At this stage, the spin temperature starts 
veering towards Tcmb until it is completely dominated by it. 
At lower redshifts the first astrophysical objects that heat and 
ionize the IGM couple T spin to the gas. Here, broadly speaking, 
there are two possible histories, one in which T sp i n couples to 
the gas as it heats up where it obtained temperature greater than 
Tcmb (red solid line). In the other possible evolution the spin 
temperature couples to the gas much before the kinetic temper- 
ature exceeds that of the CMB (red dashed line) (Thomas et al. 
2010; Baek et al. 2009). In the former case the 21 cm radiation, 
after decoupling from the CMB at z i 30, is seen only in emis- 
sion, whereas in the latter case it is seen initially in absorption 
and only at later stages in emission. 
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Fig. 2. The global evolution of the CMB (blue line), gas (green line) 
and spin (red solid line and red dashed line) temperatures as a function 
of redshift. The CMB temperature evolves steadily as 1 +z whereas the 
gas and spin temperatures evolve in a more complicated manner (see 
text for detail). 

4.2. The Brightness Temperature 

In radio astronomy, where the Rayleigh- Jeans law is usually 
applicable, the radiation intensity, /(v) is expressed in terms of 
the brightness temperature, so that 

2v 2 

/(v) = -^k B T b , (2) 

where v is the radiation frequency, c is the speed of light and 
kg is Boltzmann's constant (Rybicki, & Lightman 1986). This 
in turn can only be detected differentially as a deviation from 
Tcmb- The predicted differential brightness temperature devia- 
tion from the cosmic microwave background radiation is given 
by Field (1958, 1959); Ciardi & Madau (2003), 

5T h = 28mK (1 + 5) x m ( 1 - x 

\ * spin J 

x (o^&)^)(^)' (3) 

where h is the Hubble constant in units of 100 km s^'Mpc -1 , 5 
is the mass density contrast, xhi is the neutral fraction, and Q. m 
and Q.b are the mass and baryon densities in units of the critical 
density. Note that the three quantities, 5, xhi and T sp j n , are all 
functions of 3D position. 

One also should add to this equation the contribution of redshift 
distortion. In the linear regime of gravitational instability this 
component is simple to add (Kaiser 1987). However, in the 
quasi-linear regime this term is not well understood and should 
be explored in more detail. 

Equation 3 gives the cosmological signal we are after and 
clearly shows the complication and, at the same time, the 
wealth of information that could be extracted from the EoR 
signal. For example, according to recent simulations (e.g., Iliev 
et al. 2008; Thomas et al. 2009) the ionized fraction of the 



IGM at z ~ 1 1 is very small. If the spin temperature is much 
larger than Tcmb the brightness temperature is proportional to 
the cosmological density and can provide a clean probe of the 
primordial density field. However, at later stages one can probe 
the evolution of the ionized fraction as a function of redshift 
which provides information about the ionization sources and 
their clustering properties and so on. Figure 3 shows a typical 
distribution of the differential brightness temperature as pre- 
dicted by the recent Thomas et al. (2009) simulations. 

One should also note that in most simulations the spin temper- 
ature is assumed to be much larger than the CMB temperature, 
namely term (l - TcmbIT spin) in eq. 3 is 1. As figure 2 shows, 
this is a good assumption at the later stages of reionization, 
however, it is not valid at the early stages. Modeling this ef- 
fect is somewhat complex and requires radiative transfer codes 
that capture the Lyman-ff line formation and multifrequency 
effects, specially those coming from energetic photons (Baek 
etal. 2009; Thomas et al. 2010). 

5. The redshifted 21 cm Observation 

The current generation of radio telescopes (LOFAR b , MWA C , 
GMRT, and 2 1 CMA d ) will capture the lower redshift part of the 
5Tb evolution (z ~ 12). LOFAR, for example have two obser- 
vational bands, high band and low band. The high band array 
is expected to be sensitive enough for measuring the redshifted 
21 cm radiation coming from the neutral IGM within the red- 
shift range of z= 1 1 .4 ( 1 1 5 MHz) to z=6 (203 MHz), with a res- 
olution of 3-4 arcminutes and a typical field of view of ~ 120 
square degrees (with 5 beams) and a sensitivity of the order 
of 80 mK per resolution element and 1 MHz frequency band- 
width. At frequencies below the FM band, probed by the low 
band array, the LOFAR sensitivity drops significantly and the 
sky noise increase dramatically that detection of HI signals at 
these frequencies is beyond reach with LOFAR (Harker et al. 
2010; Jelic etal. 2008; Labropoulos et al. 2009). The main and 
most challenging task of the LOFAR EoR project is to extract 
the cosmological signal from the data and interpret it. 

In the future SKA can significantly improve on the current in- 
struments in two major ways. Firstly, it will have at least an or- 
der of magnitude higher signal-to-noise which will allow much 
better statistical detection of the EoR and give us access to 
the Universe's Dark Ages which corresponds to much higher 
redshifts (z « 35) hence providing crucial information about 
cosmology which non of the current telescope hopes to have. 
Secondly, SKA will have a resolution better by a factor of few 
at least relative to the current telescopes. These two advantages 
will not only improve on the understanding we gain with cur- 
rent telescopes but give the opportunity to address a host of 
fundamental issues that current telescope will not be able to 
address at all. Here I give a few examples: 1- Due to the lim- 
ited resolution and poor signal-to-noise, the nature of the ion- 
izing sources is expected to remain poorly constrained; 2- the 

b http://www.lofar.org & http://www.astro.rug.nl/~LofarEoR 
c http://www.mwatelescope.org/ 
d http://21cma.bao.ac.cn/ 
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Fig. 3. Position-redshift slices through the image-frequency volume of the reionization simulations of Thomas et al. (2009). 



mixing between the astrophysical effects and the cosmologi- 
cal evolution is severe during the EoR but much less so dur- 
ing the Dark Ages an epoch beyond the reach of the current 
generation of telescope but within SKA reach; 3- at redshifts 
larger than 30 the 21 cm could potentially provide very strong 
constraints, much more so than the CMB, on the primordial 
non-gaussianity of the cosmological density field which is es- 
sential in order to distinguish between theories of the very early 
Universe (e.g., distinguish between different inflationary mod- 
els). 

In section 4 we discussed the cosmological 21 cm signal and 
showed that it is expected to be of the order of x 10mK, 
However, the detectable signal in the frequency range that cor- 
responds the epoch of reionization is composed of a number 
of components each with its own physical origin and statistical 
properties. These components are : 1- the 21 cm signal coming 
from the high redshift Universe. 2- galactic and extra-galactic 
foregrounds, 3- ionospheric influences, 4- telescope response 
effects and 5- noise (see Figure 1). Obviously, the challenge of 
the upcoming experiments is to distill the cosmological signal 
out of this complicated mixture of influences. This will cru- 
cially depend on the ability to calibrate the data very accurately 
so as to correct for the ever changing ionospheric effects and 
variation of the instruments response with time. In this section 
I present the various non cosmological components and briefly 
discuss their properties. 

5.1 . Station configuration and uv coverage 

The low frequency arrays must be configured so that they have 
a very good uv-coverage. This is crucial to the calibration effort 
of the data where a filled uv plane is important for obtaining 
precise Local (Nijboer et al. 2006) and Global (Smirnov & 
Noordam 2004) Sky models (LSM/GSM; i.e. catalogues of 
the brightest, mostly compact, sources in and outside of the 
beam, i.e. local versus global). It is also crucial for the ability 
to accurately fit for the foregrounds (Harker et al. 2009b; Jelic 



et al. 2008) and to the measurement of the EoR signal power 
spectrum (Bowman et al. 2006; Harker et al. 2010; Hobson & 
Maisinger 2002; Santos et al. 2005) 

The uv coverage of an interferometric array depends on the 
layout of the stations (interferometric elements), their number 
and size as well as on the integration time, especially, when the 
number of stations is not large enough to have a good instanta- 
neous coverage. 

For a given total collecting area one can achieve a better uv 
coverage by having smaller elements (stations). For example 
LOFAR has chosen to have large stations resulting in about 
* 10 3 baselines in the core area. Such a small number of base- 
lines needs about 5-6 hours of integration time per field in or- 
der to fill the uv plane (using the Earth's rotation). In compar- 
ison, MWA which has roughly 1/2 the total collecting area of 
LOFAR chose to have smaller stations with about « 10 5 base- 
lines resulting in an almost instantaneous full uv coverage. 
The decision on which strategy to follow has to do with a num- 
ber of considerations that include the ability to store the raw 
visibilities, hence, allowing for a better calibration, an accept- 
able noise level for both the foreground extraction needs as well 
the power spectrum measurement (see the following sections; 
§ 5.2 and § 5.3). A compromise between these issues as well 
the use of the telescopes for other projects is what drives the 
final lay out of the antennas. 

5.2. Noise Issues 

The thermal noise level for a given visibility, i.e., uv point, is, 

AV(u,v)*^^L, (4) 
edASBt 

where T sys is the system temperature, e is the efficiency, dA 
is the station area, B is the bandwidth and t is the observation 
time (see e.g., Morales 2005). This expression is simple to 
understand, the more one observes, either in terms of integra- 
tion time, frequency band width or station collecting area the 
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Fig. 4. This figure shows how two different experiments might sample 
an annulus in uv. The size of uv point is given by the station (interfer- 
ometric element) size, larger station (left panel) has a larger footprint 
relative to the smaller station case (right panel) in uv plane; the foot- 
print is shown by the purple circles. Even though the sampled area in 
the two cases might be the same, the fact that smaller stations sample 
the annulus more results in an increased accuracy in their estimation 
of the power spectrum. 



less uncertainty one has. Obviously, if the signal we are after 
is well localized in either time, space or frequency the relevant 
noise calculation should take that into account. 

In order to calculate the noise in the 3D power spectrum, the 
main quantity we are after, one should remember that the fre- 
quency direction in the observed datacube is proportional to 
the redshift which in turn can be easily translated to distance 
whereas the u and v coordinates are in effect Fourier space co- 
ordinates. Therefore, to calculate the power spectrum first first 
shoud Fourier transform the data cube along the frequency di- 
rection. Following Morales (2005) I will call the new Fourier 
space coordinate rj (with dr\ resolution) which together with u 
and v define the Fourier space vector u - {«, v, rf\. From this, 
one can calculate the noise contribution to the power spectrum 
at a given \u\, 



2N7 l AT 

ZJV beam 1 v cell 



- I /2 I^BTsys \ 

\ edAdrj J 



1 



Bn(\u\)t 



(5) 



where Nb eam is the number of simultaneous beams that could be 
measured, N ce u is the number of independent Fourier samplings 
per annulus and n(\u\) is the number of baselines covering this 
annulus (Morales 2005). Note that n(\u\) is proportional to 
square number of stations, hence, n(\u\)dA 2 is proportional to 
the square of the total collecting area of the array regardless 
of the station size. This means that the noise power spectrum 
measurement does not depend only on the total collecting area, 
band width and integration time, it also depends the number of 
stations per annulus. This is easy to understand as follows, the 
power in a certain Fourier space annulus is given by the vari- 
ance of the measured visibilities in the annulus which carries 
uncertainty proportional to the inverse square root of number of 
points. This point is demonstrated in figure 4. As an example, 
MWA will have the advantage of having many samplings of the 
uv plane (100 times more than LOFAR) whereas LOFAR will 
have the advantage of simultaneous multi-beams and somewhat 
larger collecting area. 
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Fig. 5. A figure showing the various cosmological and galactic com- 
ponents that contribute to the measured signal at a given frequency. 
The slices are color coded with different tales owing to the vast differ- 
ence between the range of brightness temperature in each component, 
however the figure shows the rms of the galactic foregrounds, extra 
galactic foregrounds and cosmological signal 



5.3. The Foregrounds 



The foregrounds in the frequency regime (40 - 200MHz) are 
very bright and dominate the sky. In fact the amplitude of the 
foreground contribution, T^, at 150MHz is about 4 orders of 
magnitude larger than that of the expected signal. However, 
since we are considering radio interferometers the important 
part of the foregrounds is that of the fluctuations which reduces 
the ratio between the them and the cosmological signal to about 
2-3 orders of magnitude, which is still a formidable obstacle to 
surmount. 

The most prominent foreground is the synchrotron emission 
from relativistic electrons in the Galaxy, this source of con- 
tamination contributes about 75% of the foregrounds. Other 
sources that contribute to the foreground are radio galaxies, 
galaxy clusters, resolved supernovae remnant, free-free emis- 
sion provide 25% of the foreground contribution (see Shaver 
et al. 1999). Figure 5 shows simulated foreground contribu- 
tion at 120 MHz taken into account all the foreground sources 
mentioned. 

As many studies have shown, the very smooth structure of the 
foreground sources along the frequency direction will enable 
disentangling their contribution from that of the cosmological 
signal. The foregrounds are normally fitted by some procedure 
(e.g., polynomial fitting (Jelic et al. 2008), or more advanced 
non parametric methods (Harker et al. 2009b)) in order to re- 
cover the EoR cosmological signal; Figure 6 shows how suc- 
cessful such a recovery is. 
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Fig. 6. Statistical detection of the EoR signal from the LOFAR-EoR 
data maps that include diffuse components of the foregrounds and re- 
alistic instrumental noise (cr„ o , se (150 MHz) = 52 mK). Black dashed 
line represents standard deviation (cr) of the noise as a function of 
frequency, cyan dashed line cr of the residuals after taking out smooth 
foregrounds component and red solid line the cr of original EoR signal. 
The grey shaded surface represents the 90% of detected EoR signals 
from 1000 independent realisations of the noise, where cyan dashed 
line is mean of detected EoR signal. Note that the y-axis is in logarith- 
mic scale (Jelic et al. 2008). 
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Fig. 7. Power spectra of the cosmic signal (blue solid line), the noise 
(red dashed line), the residuals (magenta dotted line) and the extracted 
signal (black points with error bars) at three different redshifts.for the 
case when the uv coverage is frequency-dependent, we have 300 hours 
of observation per frequency channel with a single station beam, and 
the foreground fitting is done using Wp smoothing in Fourier space 
(Harkeretal. 2010). 



6. The Statistics of the observed cosmological 
signal 

6.1. The 21 cm Power Spectrum 

One of the main aims of the EoR projects is to measure the 
power spectrum of variations in the intensity of redshifted 21- 
cm radiation from the EoR. As shown in Equation 3 the power 
spectrum depends on a number of astrophysical and cosmolog- 
ical quantities. The sensitivity with which this power spectrum 
can be estimated depends on the level of thermal noise (Eq. 5) 
and sample variance, and also on the systematic errors arising 
from the extraction process, in particular from the subtraction 
of foreground contamination. In the LOFAR case, for example, 
we model the extraction process using realistic simulations of 
the cosmological signal, the foregrounds and noise. In doing 
so we estimate the sensitivity of the LOFAR EoR experiment 
to the redshifted 21 cm power spectrum. Detection of emission 
from the EoR should be possible within 360 hours of obser- 
vation with a single station beam. Integrating for longer, and 
synthesizing multiple station beams within the primary (tile) 
beam, then enables us to extract progressively more accurate 
estimates of the power at a greater range of scales and redshifts 
(see Figure 7 Harker et al. 2010). 



6.2. High order statistics: Skewness and Kurtosis 

Figure 8 showes the PDF of the brightness temperature at four 
different redshifts; the PDF is clearly nongaussian at all four 
cases. High order moments, like the skewness, as a function of 
redshift could be a useful tool for signal extraction in the pres- 
ence of realistic overall levels of foregrounds and noise. Harker 
et al. (2009a) have shown that the cosmological signal, under 
generic assumptions, have a very characteristic pattern in the 
skewness as a function of redshift (Figure 9). At sufficiently 
high redshifts the signal is controlled by the cosmological den- 
sity fluctuations which, in the linear regime, are Gaussian. At 
lower redshifts, and as nonlinearity kicks in, the signal starts 
getting a slightly positive skewness. As the ionization bubbles 
begin to show up the skewness starts veering towards until it 
crosses it to the negative side when the weight of the ionized 
bubbles becomes more important than the high density outliers 
-note high density outliers are likely to ionize first- but the dis- 
tribution is still dominated by the density fluctuations. At lower 
redshift the bubbles dominate the PDF and the neutral areas 
become the outliers giving rise to a sharp positive peak to the 
skewness. At redshift around 6 the instrument noise assumed 
to be Gaussian dominates driving the skewness again towards 
zero. Exploiting this characteristic behavior might allow us to 
pick up the cosmological signal with this high order statistic. 
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Fig. 8. The distribution of 6T\, in a certain cosmological simulation of 
reionization (Iliev et al. 2008) at four different redshifts, showing how 
the PDF evolves as reionization proceeds. Note that the y-axis scale in 
the top two panels is different from that in the bottom two panels. The 
delta-function at ST b = grows throughout this period while the rest 
of the distribution retains a similar shape. The bar for the first bin in 
the bottom-right panel has been cut off: approximately 58 per cent of 
points are in the first bin at z = 7.78 (Harker et al. 2009a). 




Fig. 9. Skewness of the fitting residuals from data cubes with un- 
correlated noise, but in which the residual image has been denoised 
by smoothing at each frequency before calculating the skewness. 
The three lines correspond to results from three different simulations 
(Thomas et al. 2009; Iliev et al. 2008). Each line has been smoothed 
with a moving average (boxcar) filter with a span of nine points. The 
grey, shaded area shows the errors, estimated using 100 realizations of 
the noise (see Harker et al. 2009a) 



7. Conclusions 

The imminent availability of observations of redshifted 21 cm 
radiation from the Universe's Dark Ages and the EoR will 
be one of the most exciting development in the study of cos- 
mology and galaxy and structure formation in recent years. 
Currently, there are a number of instruments that are designed 
to measure this radiation. In this contribution I have argued that 
despite the many difficulties that face such measurements they 
will provide a major breakthrough in our understanding of this 
crucial epoch. In particular current radio telescopes, such as 
LOFAR, will be able to provide us with the global history of 
the EoR progression, the fluctuations power spectrum during 
the EoR, etc., up to z ~ 11. These measurements will usher the 
study of the high redshift Universe into a new era which will 
bridge, at least in part, the large gap that currently exists in ob- 
servation between the very high redshift Universe (z ~ 1100) 
as probed by the CMB and low redshift Universe (z ~ 6). 

Although the current generation of telescopes have a great 
promise they will also have limitations. For example they have 
neither the resolution, the sensitivity nor the frequency cover- 
age to address many fundamental issues, like the nature of the 
first sources. Crucially, they will not provide a lot of informa- 
tion about the Dark Ages which is only accessible through very 
low frequencies in the range of 40 - 120 (z ~ 35 - 1 1). 

Fortunately, in the future SKA can improve dramatically on the 
current instruments in three major ways. Firstly, it will have at 
least an order of magnitude higher signal-to-noise which will 
allow much better statistical detection of the EoR. Secondly, 
it will give us access to the Universe's Dark Ages which cor- 
responds to much higher redshifts (z ~ 30) hence providing 
a truly pristine probe of cosmology. Thirdly, SKA will have, 
at least, a factor of few better resolution in comparison with 
the current telescopes, thus better constraining the dominant 
ionization sources. In summary these three advantages will not 
only improve on the understanding we gain with current tele- 
scopes but give the opportunity to address a host of fundamen- 
tal issues that current telescope will not be able to address. 

The next decade will be extremely exciting for studying the 
high redshift Universe, especially as these radio telescope grad- 
ually come online, starting with LOFAR and MWA. As they 
promise to resolve many of the puzzles we have today pertain- 
ing to the formation and evolution of the first object, cosmol- 
ogy and the physical process in the high redshift intergalactic 
medium. 

Acknowledgements. Many of the results shown here have been ob- 
tained in collaboration with the members of the LOFAR EoR project 
whose contribution I would like to acknowledge. 
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